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Abstract Applying the Horsky-Mitskievitch conjecture to the empty space solutions 
£\J , of Morgan and Morgan due to the gravitational field of a finite disk, we have obtained 

the corresponding solutions of the Einstein-Maxwell equations. The resulting expres- 
| sions are simply written in terms of oblate spheroidal coordinates and the solutions 

represent fields due to magnetized static thin disk of finite extension. Now, although the 
solutions are not asymptotically flat, the masses of the disks are finite and the energy- 
momentum tensor agrees with the energy conditions. Furthermore, the magnetic field 
and the circular velocity show an acceptable physical behavior. 
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1 Introduction 

The study of axially symmetric solutions of the Einstein and Einstein-Maxwell field 



(N 

equations corresponding to disklike configurations of matter, apart from its purely 
\^») | mathematical interest, has a clear astrophysical relevance. Indeed, thin disks can be 

used to model accretion disks, galaxies in thermodynamical equilibrium and the su- 
perposition of a black hole and a galaxy. Disk sources for stationary axially symmetric 
spacetimes with magnetic fields are also of astrophysical importance mainly in the 
study of neutron stars, white dwarfs and galaxy formation. 

Exact solutions that have relativistic static thin disks as their sources were first 
studied by Bonnor and Sackfield [T] and Morgan and Morgan [2j[3] . Subsequently, 
several classes of exact solutions corresponding to static 4,5 6,7,8,9, 10 1112 13 and 
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stationary 14,15,16,17 thin disks have been obtained by different authors, and the su- 
perposition of a static or stationary thin disk with a black hole has been considered [18] 
[Tgil20ll21ll22ll23ll24ll25ll26| . Relativistic disks embedded in an expanding FRW universe 
have been studied in [27], perfect fluid disks with halos in [2S], and the stability of 
thin disks models has been investigated using a first order perturbation of the energy- 
momentum tensor in [29]. On the other hand, thin disks have been discussed as sources 
for Kerr- Newman fields [301131] . magnetostatic axisymmetric fields [35], and confor- 
mastatic and conformastationary metrics 33,34,35 . Also, models of electrovacuum 
static counterrotating dust disks were presented in [36] , charged perfect fluid disks were 
studied in [37] , and charged perfect fluid disks as sources of static and Taub-NUT-type 
spacetimes in [3"8"ll3"5] . 

Now, the thin disks with magnetic fields presented at references 30,31,32 were 
obtained by means of the well known 'displace, cut and reflect' method in order to 
introduce a discontinuity at the first derivative of one otherwise smooth solution. The 
result is a solution with a singularity of the delta function type in all the z = hypersur- 
face and so can be interpreted as an infinite thi disk. On the other hand, solutions that 
can be interpreted as thin disks of finite extension can be obtained if a proper coordi- 
nate system is introduced. A coordinate system that adapts naturally to a finite source 
and presents the required discontinuous behavior is given by the oblate spheroidal co- 
ordinates. Some examples of finite thin disks from vacuum solutions expressed in these 
coordinates can be found in references [1|'2HI7|, and from electrovacuum solutions in 
reference [35] . 

According to the above considerations, the purpose of our paper is to present a new 
infinite family of exact solutions of the Einstein-Maxwell equations for axially sym- 
metric spacetimes. The solutions are obtained by the use of the Horsky-Mitskievitch 
conjecture [40] . which prescribes a quite close connection between isometries of vac- 
uum spacetimes (seed metrics) and the electromagnetic four-potential of a generated 
Einstein-Maxwell fields. We take the Morgan and Morgan metric disk [2] as the seed 
solution. The generated solutions describe a family of magnetized finite thin disks, 
which is the magnetized version of the family of relativistic static Morgan and Morgan 
disks. 

The plan of our paper is as follows. First, in Section]!?] the Weyl-Lewis-Papapetrou 
line element is considered and Einstein-Maxwell equations in cylindrical coordinates 
are introduced. The procedure to obtain magnetovacuum static, axially symmetric rel- 
ativistic thin disks without radial pressure is also summarized in this section. Section 
[3] introduces the Horsky-Mitskievitch generating conjecture used to obtain a general 
kind of Weyl-Lewis-Papapetrou spacetime with magnetic field. In Section [4] we put 
forward a solution of the Einstein-Maxwell equations describing an infinite family of 
finite static magnetized thin disks by the use of the Horsky-Mitskievitch (HM) conjec- 
ture. For this purpose, the well known Morgan and Morgan metric disk [2J is employed 
as seed solution & s . The analysis of the physical behavior of solutions is presented in 
Section [5] where the asymptotic behavior of the solutions is examined. Then, we study 
the behavior of the corresponding energy density, current density, circular velocity and 
azimuthal pressure. Finally, in Section [6] we conclude with a discussion of our results. 
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2 Einstein-Maxwell Equations and Thin Disks 

The vacuum Einstein-Maxwell equations, in geometrized units such that c = SnG = 
HO = eo — 1) can be written as 

G a b = T ab , (la) 
F ab . b = 0, (lb) 

with the electromagnetic energy-momentum tensor given by 

T a b = FacF b C — -g a bFcd,F Cd , (2) 

where 

F a b = A b ,a — A a .b 

is the electromagnetic field tensor and A a is the electromagnetic four potential. The 
line element of a static vacuum spacetime can be written in the Weyl-Lewis-Papapetrou 
(WLP) form gl] 

ds = — e dt + e [r dip +e (dr + dz )J, (3) 

where t,ip,r,z are the usual cylindrical coordinates: — oo < t,z < oo, < r, < 
ip < 2-k, the metric functions <& and /l depend only on r and z. 

The solutions of the Einstein-Maxwell equations corresponding to a disklike source 
are even functions of the z coordinate. Therefore, they are continuous functions every- 
where but their first z-derivatives discontinuous at the disk surface. Consequently, in 
order to obtain the energy-momentum tensor and the current density of the source, 
the jump across the disk of the first z- derivatives of the metric tensor is expressed as 

Kb = [9ab,z] = 2 9ab, z \ z=0+ , (4) 
and the jump across the disk of the electromagnetic field tensor is expressed as 

[F za ] = [Aa, z ]=2A a , z \ z=0+ , (5) 

where the reflection symmetry of the functions with respect to z = has been used. 

Then, using the distributional approach 42.431H1] or the junction conditions on 
the extrinsic curvature of thin shells 45,46,47 , the Einstein-Maxwell equations yield 
an energy-momentum tensor as 

T ab = T ab e ^ + T ab^ _ + Q^g^ (g) 

and the current density as 

J a = I a S(z), (7) 

where 6(z) and 8(z) are, respectively, the Heaviside and Dirac distributions with sup- 
port on z = 0. Here T± h are the electromagnetic energy-momentum tensors as they are 
defined by © in the z > and z < regions, respectively, whereas 



-.a fja Z £ Z i ZZ ca . a Z j Z ZZ ia . iC/ ZZ c-a a2 r2\-i / \ 

h = ni b s b~ b k+9 h~9 h+b c {g h~9 ( 8 ) 



gives the part of the energy-momentum tensor corresponding to the disk source, and 

I a = [F az ] (9) 
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is the contribution of the disk source to the current density. The surface energy- 
momentum tensor, S a i,, and the surface current density, j a , of the disk, can be obtained 
through the relations 

}ab S(z) ds n = e A ~® Q ab , (10) 



and 

ja = 



J I a 6(g) ds n = e A -*I a , (11) 



where ds n ~ ^Jgzz dz is the physical measurement of length in the normal direction of 
the disk. 

Now, we choose the magnetic potential as 

A a = (0,A(r,z), 0,0), (12) 
so that we have a pure axially symmetric magnetic field 

F^pz = ^,2j Frtp = -A,r- (1^) 
Thus, for the metric ([3]), the only non-zero components of S% and j a are 

S§= 2e®~ A {A, z - 2$, z } , 

S{ = 2e 4 '- A A,z , (14) 
jtp = ~2e 

where all the quantities are evaluated at z = + . 
With the orthonormal tetrad 

e (6 a ) = {V a ,W a ,X a ,Y a }, (15) 

where 

V = e~*(l, 0,0,0), (16a) 

W a = —(0,1,0,0), (16b) 
r 

X tt = e*- A (0, 0,1,0), (16c) 

y a = e*~ A (0, 0,0,1), (16d) 

we can write the surface energy-momentum tensor of the disk in the canonical form as 

S ab = eV a V b +pW a W b . (17) 

Here e and p are the energy density and the azimuthal pressure of the disk, respectively. 
We also consider the mass density on the surface of the disk, defined as 

V = e + p. (18) 

In the same way, the current density on the disk can be written as 

j a =jW a . (19) 
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As we can see the charge density on the disk surface is identically equal to zero. This 
can be checked by using Eqs. (1111) . (1121) and (1191) . The electromagnetic field is of a 
magnetic type, which can be demonstrated by the electromagnetic invariant 

F *bF = fl 7 2 [A* r + A%] > 0. (20) 

We consider now, on the basis of Refs. P2] and [35] , the possibility that the energy- 
momentum tensor S ab and the current density j a can be interpreted as the superpo- 
sition of two counterrotating fluids. In order to do this, we cast 

s ab = e+U lu b + + e _U a LU h _, (21a) 



f = a+U+ +a-U°l, (21b) 

where 

e+ = e_ = (e - p)/2, (22a) 



=^v?- 1, (22b) 

are the energy densities and charge densities of the two counterrotating fluids. The 
counterrotating velocity vectors are given by 

V a ± UW a 



where 



U 2 = H < 1 (24) 



is the counterrotating tangential velocity. Therefore, we have two counterrotating 
charged fluids with equal energy densities and equal but opposite charge densities. 



3 The Horsky-Mitskievitch conjecture and the WLP metric 

3.1 The seed WLP metric 

The WLP line element of a static vacuum spacetime can be written as [41] 

j 2 2<5 S ,,2 . -2* s r 2 , 2 . 2A S / , 2 , , 2\i , n ri 

as — — e at + e [r dip + e (ar + az )\, (25) 

where {t, tp, r, z} are the usual cylindrical coordinates: — oo < t, z < oo, < r, < ip < 
2tt. The metric functions $ s and A s depend on r and z. This metric admits two Killing 
vectors: 



it = (1,0,0,0), 

& = (0,1,0,0), 



(26a) 
(26b) 
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which, in the vectorial space di are, 

€i = (27a) 

^ = ( 27b ) 

or, in terms of their components are 

t{ = 4, (28a) 

^ = 4. ( 28b ) 

where 8% is the usual Kronecker delta tensor. 
In terms of the one-forms dt and d(p we have 

£t = -e 2 * s dt, (29a) 

£ v = r 2 e~ 20 'dAp, (29b) 

where all the orthonormal bases employed are always chosen as a generalization of the 
set 

w<°> = e* s dt, (30a) 



^(l) _ re 



-'Ik 



dip, (30b) 
w 0O = re ^«-*.d r , (30c) 
w (s) = re ^-*-d«, (30d) 



the simplest tetrad that can be used for the Weyl solution. 



3.2 The WLP solution with a magnetic field 

The present section is devoted to the derivation of solutions of the Einstein-Maxwell 
equations via HM conjecture. Following the method used in [3S], we modify the line 
element (|25[) to the form 

— 2#, 2 

j 2 p / \ 2 1<S> S j.2 , e r , 2 . n2 2(yl»— 2 , , 2\ / Q1 \ 

as = -f(r,z) e s dt + — ^-dtp +f(r,z) e B '(dr + dz ), (31) 

f{r,zy 

that admits two Killing vectors 

& = -f(r,z) 2 e 2 *°dt, (32a) 

r 2 e -2<S 3 

£ v = -77 ^"dy, (32b) 

with /(r, z) an arbitrary function. 

The Killing vectors £ and the electromagnetic four-potential A satisfy 

*d*d£ = 0, *d*dA = 0, (33) 

where the * is the usual Hodge (star) operation. Then, the Killing vector £ v induces 
the four-potential 

A = qf(r,z)$p, (34) 
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It can be verified through standard calculations that the sourceless Einstein-Maxwell 
equations are fulfilled if 

/(r,«) = l + ci/ 1 (r,*), (36) 
where the function f\ must be a solution of the differential equation 

G ia \ a) = —R = 0, (37) 

given the fact, for electrovacuum spacetimes, the Einstein tensor, C?( a )((,), is traceless. 
As in [48] we can choose c\ = q 2 and /i = ?- 2 e~ 2 * s , and then we have the four-potential 

A = 2 2<S> d ^ (38) 

qZ r z _|_ e zv s 

the q parameter characterizing the strenght of the electromagnetic field. Under the 

conditions described above the metric set by ()31|) is reduced to 

■ 2 2* ,,2 , -2<f r 2 , 2 , IA/j 2 , , 2\ n / rv\ 

as = — e at + e |r ac/5 + e (or + az )J, (39) 

where, /e* s = e # and /V s = e A . 

Using (|14|) . (|15|) and (| 1 Tp we find out the surface energy density 

4g 2 (r) e 3<l ' 3 ~ yl3 (l- r<P s , r )<P 
(e 2 * + q 2 r 

the azimuthal pressure 

_ 4r£(r)e 3 * s - As [2g 2 r- + .E(r)<£ s , r ]i> 



'_2<2> , „2 T .2>3 ' V ; 



the surface mass density 



the surface current density 



and the circular velocity 



( e 2* fl +? 2 r 2)3 



4S(r)e 3 * s - yl ^ s 



(41) 



M (e 2 ^ +q 2 r 2 ) 2 ' ' (42) 



. 4gre 4 *--^g a ,s 
(e 2 * S + g 2 r 2 ) 2 ' 



U 2 = [/? 



1 + . ^ 



E(r)<P s 



(44) 



where -E'(r) = e 2 * s — a 2 ?" 2 , and 



is the circular velocity of the unmagnetized source, i.e. q = 0. 

As we can see from (|29ap . (|29b|) . ()32ap and (|32bf) . the method is a solution gen- 
erating technique that changes the norm of the two Killing vectors without adding 
twist. Besides, the possible connection of the HM conjecture with inner symmetries 
of Einstein-Maxwell equations would connect the conjecture with a set of generating 
methods elaborated by Ernst and other authors for axisymmetric fields |49] , 
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4 Magnetized Morgan-Morgan disks 

Let's restrict the previous general model to obtain a solution of the Einstein-Maxwell 
equations describing an infinite family of finite static magnetized thin disks. For this 
purpose, we use, as a seed solution (<£ s , A s ), the well known Morgan and Morgan 
metric disk (</>, A) [2], 

oo 

<t>(x,V) = ~^2 C 2n q2n{x)P2n(y); 

71=0 

where Cm are arbitrary constants that must be properly specified so that a partic- 
ular solution may be set. P n {y) are the usual Legendre polynomials and q n {x) = 
i n+1 Q n (ix). Qn(z) being the Legendre functions of second kind (see |50] and, for the 
Legendre functions of imaginary argument, [51] pag. 1328). The x and y are the oblate 
spheroidal coordinates related with the cylindrical coordinates by the relations [51] 

r 2 = a 2 {l + x 2 )(l-y 2 ), (46a) 
z = axy, (46b) 

where < x < oo and — 1 < y < 1. The disk has coordinates x = 0, < y 2 < 1 and, 
when the disk is crossed, the sign of y changes, but not its absolute value. 

We use the constants Cm as they have been determined by Gonzalez and Reina in 



mG 

Cm — — 

la 



7r 1/2 (4n + l)(2Z + l)! 



22' (2n + !)(/- n)\r{l + n+ §)g 2 n+l(0) 



for n < I and Cm = for n > I. Consequently, these solutions correspond to the 
magnetized version of the well known Morgan-Morgan disk [2], For instance, the first 
three members of the family of the seed metric functions are given by 

<h.(x,y) = -^[cot" 1 x + A(3y 2 - 1)], (47a) 

2 (x,y) = -^[cot- 1 a; +^(3y 2 -l)+ B(35y 4 - 30y 2 + 3)], (47b) 

<j> z {x, y) = - — [cot' 1 x + ^(3y 2 - 1) + (47c) 
a b 

^(35y 4 - 30y 2 + 3) + C(23l/ - 315y 4 + W5y 2 - 5)], 



where 



A = ^-[(3a; 2 + 1) cot 1 x - 3x], 

B = -^r[(35a; 4 + 30x 2 + 3) cot" 1 x - 35x 3 - %-x\. 
448 3 

C = — — \(231x 6 + 315a; 4 + 105a; 2 + 5) cot -1 x - 231a; 5 - 238a; 
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-x\ 



8448 LV ' '5 

which similar, but more involved, expressions for greater values of I. 

On the other hand, according to H40[) . the energy density of the magnetized disk 
can be written as 

, _ 4E 2 e 3 ^^[y + (l-y 2 ) ( j >n ^ n , x 

ay 2 [e 2 ^ +q 2 a 2 (l-y 2 )] 3 ' 1 1 
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with 



while, according to (|41|l and (|42|l the pressure, p, and mass density, fx, in the surface 
of the disk are 



Pn 



4E(1 - j/V^-^.spgVy - <j>n,yE] 
ay 2 {e 2 <l>n +q 2 a 2 (l-y 2 )] 3 



(50) 



and 



Hn = 



(51) 



respectively, and, according to (|43p and (|44f) . we have the following surface current 
density and circular velocity 



4 g (i- y Y/V^-^ 

y[e 2c t>™ + q 2 a 2 (l - y 2 )] 2 ' 



(52) 



and 



ul = t/ s 2 



2g 2 a 2 j/ 



E(j> n ,y 



(53) 



respectively, where 



y + (i-y 2 )<f>n,y' 



is the circular velocity of the unmagnetized source, i.e. q — 0. The non zero components 
of the magnetic field on the surface of the disk, B z — -~A VtT and B r = A v ^ z , are 



B z = 



2qa(l-y 2 ) 1 / 2 e 2 ^[y+(l-y 2 )<j ) r h y] 
y\e 2 <t>™ + q 2 a 2 (l - y 2 )} 2 



(54) 



and 



Br — 



2ga(l - y 2 )e 2< ^», a 

'y[ e 2<j,n +g 2 a 2( 1 _ J/ 2)]2' 



(55) 



where we have used the magnetic potential of ()35[1 , and all the quantities are evaluated 
on the disk. 
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5 The behavior of the solutions 

In order to elucidate the behavior of the different particular models, firstly we introduce 
the dimensionless energy surface density on the disks, defined for < f < 1 as e n {r) — 
ae n (r), where the dimensionless radial variable r = af has been introduced. In figure[T] 
the dimensionless surface energy densities l n for the models corresponding to n = 1, 2 
and 3 are depicted. In each (f) for < r < 1 with m — 0.10 for different 

values of the parameter q — aq are shown. First, we set q = 0.30, the bottom curve 
in each graphic, and then 0.20, 0.10 and q — (dotted curve), the vacuum case. In 
all the cases, it can be seen that the energy density is everywhere positive and that 
it vanishes at the edge of the disk. The disks with higher values of n and q show an 
energy distribution that is more concentrated in the center and lower at the edge. It 
can be observed that the energy density in the central region of the disk decreases with 
the presence of a magnetic field. 

The figure [5] shows the dimensionless azimuthal pressure, p n (r) = ap n (r), for the 
models corresponding to n = 1, 2 and 3. As we see, in all the cases, the pressure is 
everywhere positive and rapidly increases as we move away from the center of the disk. 
It reaches a maximum and, later, it rapidly decreases. It can be observed that the 
presence of a magnetic field increases the pressure everywhere on the disk. Otherwise, 
the maximum increases and moves away from the disk edge as n increases. It can be 
also observed that the surface current density, represented in figure [4j has a behavior 
similar to that of the pressure. 

In order to illustrate graphically the behavior of the circular velocities or rotation 
curves, we introduce the dimensionless quantity U n (r) = aU n (r). In figure [3] we plot 
the dimensionless rotation curves for the models corresponding to n = 1, 2 and 3. 
The circular velocity corresponding to n — 1 is a monotonously increasing function of 
the radius. On the other hand, for n — 2 and n — 3, the circular velocity increases 
from a value of zero at the center of the disks until it attains a maximum at a critical 
radius and then, it decreases to acquire a finite value at the edge of the disk. It can 
also be seen that the value of the critical radius increases as the values of n decreases 
and q increases. From which, we can conclude that the magnetic field decreases the 
circular velocity everywhere on the disk. We study these solutions with other values of 
the parameters as well, but, in all the cases, we find a similar behavior. We see that 
the magnitude of m and q are limited uniquely by the condition that the magnitude 
of the velocities of the disk cannot exceed the speed of light. The explicit calculation 
of the Kretschmann invariant R a bcdR abcd is too lengthy to be reproduced here. A 
direct computation by the use of REDUCE shows that the first member of the family 
of the magnetized Morgan-Morgan disk, n = 1, has a singularity at the rim of the 
disk where the Kretschmann invariant becomes infinite although its mass density is 
finite everywhere, when all the n > 1 disks are regular. Obviously, this property in the 
curvature is, of course, inherited from the seed Morgan and Morgan metric disk (see 
[53]). 

The solution of ()39|) in terms of the seed solution (Morgan and Morgan metric disk, 
&s — <f> ) can t> e written in the form 

<£ = ln[e 2 * s + q 2 r 2 ] -<P S . (56) 

Therefore, when we take r — > oo clearly we see that, unlike the Morgan and Morgan 
seed solution, generated spacetimes are not asymptotically flat. Consequently, although 
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the technique described here produces a new class of solution of Einstein-Maxwell's 
equations and all the physical quantities exhibit an excellent behavior, the method 
offers very little insight into the physical implications of the metrics. The reason of 
this non-flatness far away from the disk should be that the magnetic field itself does 
not decay at infinity as it appears to be easily expressible from (|20|l . On the other 
hand, this still appears to be a viable solution of the Einstein-Maxwell system but, 
apparently, corresponds to magnetized disks immersed in a Melvin magnetic universe. 



6 Concluding remarks 

We have presented an infinite family of new exact solutions of the vacuum Einstein- 
Maxwell equations for static and axially symmetric spacetimes. The solutions describe 
an infinite family of magnetized finite thin disks, the magnetized version of the family 
of Morgan and Morgan relativistic thin disks [2j. The first member of the family of 
the magnetized Morgan-Morgan disk, n = 1, has a singularity at the rim of the disk 
where the Kretschmann invariant becomes infinite, although its mass density is finite 
everywhere. Whereas all the n > 1 disks are regular. Obviously, this property in the 
curvature is, of course, inherited from the seed Morgan and Morgan metric disk. 

Unlike the Morgan and Morgan seed solution, the generated spacetimes are not 
asymptotically flat. Consequently, although the technique described here produces a 
new class of solution of Einstein-Maxwell's equations and all the physical quantities 
exhibit an excellent behavior and the energy-momentum tensor is in fully agreement 
with all the energy conditions, the method offers very little insight into the physical 
implications of the metrics. However, these are the first fully integrated exact solu- 
tions for this kind of magnetized thin disk sources. They are important in so far as 
they represent a new family of exact solutions of the Einstein-Maxwell vacuum equa- 
tions. Moreover, the outlined method may serve as a guideline to find more physically 
acceptable solutions in future works. 



Appendix 

We pass on to the definition of the Hodge * (star) operation. Let's us M a m— dimensional 
manifold endowed with a metric g. The Hodge * (star) operation is a map: f2 r (M) — > 
S7 , ~ m (M) whose action, on the basis of the vector Q r {M), is defined by |54] 

*(dx fl1 A dx^ 2 A ... A dx^) = V\9~\ »i»2.~»r dx »r+i A _ _ _ A dx » m _ 

(m — r)! rr-t-x—rm 

For a r— form 

1 



A = ■^A fllMa ... Atr dsc /l1 A dx^ 2 A ... A dx^ r G O r 



we have 



*a = cfe/^i a ... a dx^ e n r > 

r!(r — my. * + ^ 
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n=1 




0.2 0.4 0.6 0.8 1 

n=2 




0.2 0.4 0.6 0.8 1 



Fig. 1 Dimensionless energy density e n as a function of f for the first three disk models with 
n = 1, 2, 3. In each case, we plot e n (f) for < f < 1 with m = 0.10 for different values of the 
parameter q. First, we take q = 0.30, the bottom curve in each plot, and then 0.20, 0.10 and 
q = (curve with dots), the vacuum case. 
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0.14 




Fig. 2 Dimensionlcss azimuthal pressure p n as a function of f for the first three disk models 
with n = 1, 2, 3. In each case, we plot p n (f) for < f < 1 with m = 0.10 for different values of 
the parameter q. First, we take q = 0.30, the top curve in each plot, and then 0.20, 0.10 and 
q = (curve with dots), the vacuum case. 
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n=1 




Fig. 3 Dimensionless circular velocity U n as a function of f for the three first disk models 
with n = 1, 2, 3. In each case, we plot U n (r) for < f < 1 with m = 0.10 for different values 
of the parameter q. We first take q = 0.30, the top curve in each plot, and then 0.20, 0.10 and 
q = (curve with dots), the vacuum case. 
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0.18 




Fig. 4 Dimensionlcss surface current density j„ as a function of f for the first three disk 
models with n = 1,2, 3. In each case, we plot j n (f) for < f < 1 with m = 0.10 for different 
values of the parameter q. First, we take q = 0.30, the top curve in each plot, and then 0.20 
and 0.10 
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where the totally anti-symetric tensor e is 



+ 1 if (H1H2 •••Mm) 

is an even permutation of (12 . . . m), 

-1 if 0i^2 •••/Urn) 
is an odd permutation of (12 . . . m), 

otherwise, 



and "A" is the usual exterior product or wedge product. 
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